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Abstract 

We prove that the variance of the current across a characteristic is of or- 
der t 2 ^ in a stationary constant rate totally asymmetric zero range process, 
and that the diffusivity has order i 1 / 3 . This is a step towards proving uni- 
versality of this scaling behavior in the class of one-dimensional interacting 
systems with one conserved quantity and concave hydrodynamic flux. The 
proof proceeds via couplings to show the corresponding moment bounds for 
a second class particle. We build on the methods developed in [4] for simple 
exclusion. However, some modifications were needed to handle the larger 
state space. Our results translate into t 2//3 -order of variance of the tagged 
particle on the characteristics of totally asymmetric simple exclusion. 
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current fluctuations, second class particle 

2000 Mathematics Subject Classification: 60K35, 82C22 

1 Introduction 

The constant rate totally asymmetric zero range process (TAZRP) is a Mar- 
kov process that describes the motion of particles in the one dimensional integer 
lattice Z. Namely, if any particles are present at a site, then one of them jumps 
one site to the right with rate 1 independently of particles at other sites. Later 
we will give a construction of TAZRP in terms of Poisson clocks. This process is 
among the interacting particle systems introduced in [14]. Our paper examines 
the net particle current seen by an observer moving at the characteristic speed 
of the process. The characteristic speed is the speed at which perturbations 
travel in the system and can be determined e.g. via the hydrodynamic limit. 
The process is assumed to be in one of its extremal stationary distributions, 
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i.e. with independent geometrically distributed occupation variables at each site 
with parameter — j^, implying density g > of particles. We prove that the net 
particle current across the characteristic, that counts the number of particles that 
pass the observer from left to right minus the number that pass from right to left 
during time interval (0, t], has variance of order < 2 / 3 . 

The seminal papers of Baik, Deift and Johansson [1] and Johansson [9] gave 
the first rigorous proofs of such fluctuations in the last-passage version of Ham- 
mersley's process and in the totally asymmetric simple exclusion process (TASEP) 
started from jam-type initial configuration. The correct order was verified to be 
i 1 / 3 , and the limiting fluctuations were found to obey Tracy- Widom distributions 
from random matrix theory. Later a last-passage representation was also found 
for a stationary TASEP [11], and then the Tracy- Widom limit proved for the 
current across the characteristic in that setting [8]. 

Cator and Groeneboom [6] used a probabilistic approach to prove the cor- 
rect order of these fluctuations for Hammersley's process, Balazs, Cator and 
Seppalainen [3] obtained the same result for the last passage representation of 
the stationary TASEP. Then Balazs and Seppalainen in [4] proved that the cur- 
rent variance is of order t 2 / 3 in the (non-totally) asymmetric simple exclusion 
process (ASEP). 

There are generalizations of the result for the ASEP for those processes in 
which particles can jump more than one site, most recently by Quastel and Valko 
[13, 12]. 

Physical reasoning, e.g. [15], suggested universality in the sense that similar 
scaling of the current across the characteristic should occur in other systems as 
well with one conserved quantity. 

To our knowledge, the present article is the first one that generalizes the result 
of [4] towards another system in which more than one particle per site is allowed. 
Due to the connection between the TASEP tagged particle and the TAZRP par- 
ticle current, our results also apply to tagged particle fluctuations of the TASEP. 
As in [4] , our arguments are entirely probabilistic and utilize couplings of several 
processes and bounds on second class particles. Informally speaking, second class 
particles are perturbations in the system that do not disturb the motion of the 
regular particles but are influenced by the ambient system. Precise definition will 
be given after the construction of the coupled processes in Section 2. As the state 
space is larger than in the ASEP, some delicate modifications are required when 
we couple several processes together. 

Fluctuation results for asymmetric exclusion processes have also been stated 
in terms of a quantity called the diffusivity D(t). The link between current 
variance and diffusivity can be described also using the concept of the second 
class particle: the variance of the current is the expected absolute deviation of 
the second class particle, while tD(t) is the variance of the second class particle. 
For the TAZRP we also obtain the correct order i 1 / 3 for the diffusivity. 

When the observer's speed V is different from the characteristic speed, the 
fluctuations are normal, i.e. the current variance is of order t. This was proved 
by Ferrari and Fontes in [7] for the ASEP and by Balazs for some other processes 
in [2]. As a consequence we also reproduce these results concerning the normal 
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fluctuations for TAZRP. 

This paper is the continuation of [4] and might give some ideas that could lead 
to the treatment of some other models. The method seems to be robust enough 
to be able to generalize it to other, even more complicated processes. Proceeding 
in this direction towards universality is subject to future work. 

2 The constant-rate zero range process and the 
results 

Construction of the process and the second class particles 

The constant-rate totally asymmetric zero range process (TAZRP) is a Markov 
process on the state space fl = {0, 1, 2, . . . } z = Z +Z . Given a state lo = {uJi}iez € 
ft, the following jumps can happen independently at different sites: 

(2.1) (uj 1 , uj 1+1 ) — > (uji - 1, ui i+1 + 1) with rate I( Wi>0 )- 

We interpret the process as representing unlabeled particles that execute 
nearest-neighbor random walks on Z, in such a way that at a site where at least 
1 particle is present only the topmost one is allowed to jump to the right with 
rate 1. The value Ui(t) = means that site i is vacant at time t. The state of 
the entire process at time t is then lu (t) . 

A rigorous construction of this process is done by giving each site i a rate 1 
Poisson process Ni^ i+ i on the time line [0, oo). The processes {N^i + i,i e Z} 
are mutually independent, and also independent of the initial configuration w(0). 
The rule of evolution is that when N^ i+1 jumps, the topmost particle is moved 
from i to i + 1 if i is occupied by at least one particle. Thus the rates (2.1) are 
realized. 

Let /ig denote the geometric measure with mean g, i.e. n e {k} = ^pr 
on the set Z + , and let /i = /i® 2 be the i.i.d. geometric product measure with 
marginals /i g on Q. It is known that the measures {fj,^ : < g} are the extreme 
points of the convex set of invariant distributions for the process that are also 
invariant under spatial translations. 

It is convenient to embed the zero range process in a height process that 
represents a wall of adjacent columns of bricks. On top of each interval [i, i + 1] 
sits a column of bricks with height hi G Z. The entire height configuration is 
h = {hi}i € z, restricted to satisfy 

(2.2) < — hi for each i 

so that the wall slopes downward to the right. Let the Poisson processes govern 
the evolution of the heights: when A^j + i jumps add a brick on top of the column 
on + 1]. But suppress every step that leads to a violation of (2.2). 
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Given an initial particle configuration w(0), define an initial height configura- 
tion by 

(£?= i+ i"i(0) for*<0, 
hi(0)=lo fori = 0, 

l-Ej=iWi(0) for7>0. 
Let the heights evolve, and define 

u)i{t) = hi-i(t) - hi(t). 

Then this process u{t) is exactly the TAZRP constructed earlier, and the height 
increment hi(t) — /ij(0) is the net particle current across the bond (i, i + 1). 

Below we run several processes started from different initial configurations to- 
gether governed by the basic coupling in which we use common Poisson clocks for 
them. The first observation is that this coupling preserves monotonicity among 
both particle and height configurations. Ordering is defined sitewise: for parti- 
cle configurations 77 < uj_ means that r\i < u>i for each i € Z, and similarly for 
height configurations g < h if gi < hi for each i G Z. The basic coupling has the 
following property, called attractivity: 

t?(0) < w(0) => ?70) < u(t) and g(Q) < h(0) => g(t) < h(t) 

for all t > 0. 

We use the following terminology: if we have two coupled zero range processes 
•q(t) < Lu(t), then there are u>i — r\i pieces of lu — r\ second class particles at each 
site i. The joint process (r](-), w(-)) can be constructed from a two-class process: 
(i) The first class particles 77 obey the TAZRP dynamics as described earlier, (ii) 
The second class particles also obey the Poisson clocks when they can, but they 
are only allowed to jump if there are no more first class particles at site i, i.e. their 
jump rate at site i is 1(^=0) ■ I(wi>i) an d only the topmost second class particle 
is allowed to jump. We emphasize here a different behavior of the zero range 
second class particles and the ones in the simple exclusion process, i.e. a jump of 
a first class particle to site i in TAZRP does not affect the number of second class 
particles at site i. In the ASEP case, where at most one particle is allowed at one 
site, they interchange sites, implying a left jump of the second class particle even 
in the totally asymmetric case. Hence the TAZRP model second class particles 
are allowed to jump only to the right. 

Let £j € denote a configuration that has only a single particle at site i. 
For 77 <E we can legitimately define 77" 1 " = r] + 6$. In this situation we say that 
there is a single second class particle between 77" 1 " and 77 at site 0. Since the basic 
coupling conserves it, there is always a site Q{t) such that 

r L +{t) = r L (t)+5 Q{t) . 

Q(t) is the position of the second class particle at time t, which performs a nearest 
neighbor walk, influenced by the ambient process 77(-). 

It is convenient to also have the notion of a second class a??i7particle at po- 
sition Q a (t) in a process This means that Q a (t) is the location of the 
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single discrepancy between two processes LO_(t) and uT (t) that are started so that 
ujT(0) = lo(0) — 8 1 > where i = Q a (0). In the proofs we will couple more than 
two processes and this flexibility will be convenient. 

Current fluctuations and diffusivity 

Let [x] denote the first integer from x towards the origin, in other words [x] = \x\ 
(floor) when x > and [x] — \x~\ (ceiling) when x < 0. For a speed value VeR 
define 

(2.3) J (Y) (t) = h m (t), 

the height of the column over interval [ [Vt] , [Vt] + 1] at time t. The normalization 
/io(0) = implies that j( v \t) is the total net particle current seen by an observer 
moving at speed V during time interval [0,t]. One can compute the density g 
stationary expectation 

(2.4) W ji.v ){t) = _Q_ t _ e[Vt] 

by writing a martingale for ho(t) and then adding in h[ V t] (*) — h (t) which counts 
particles between sites and [Vt]. 

Our results are based on an interplay between currents and second class par- 
ticles. One key fact is the coming connection. To simplify notation we introduce 
a new measure 

^ oo 

^ V ' y = k+l 

It can be checked easily that fi is again a probability measure for any probability 
measure \i. It arises from [5, Theorem 2.2]. An easy computation gives that 

in the case of TAZRP jl e (k) = (k + (lTe) ' wmcn i s distribution 

of the sum of two independent geometrically distributed random variables with 

mean g. We will use the configuration of a coupled pair of processes started in a 

product measure where all sites have independent initial marginals [i e except the 

origin where the initial distributions are fi e and £t e + l respectively, independently 

of other sites. Denote this coupling measure of the two states w _ (0) and lo(0) by 

A • 
— e 

From now on E e and Var e denotes expectation and variance in the stationary 
process with density g (process started from /x® z ) while E and Var denotes 
expectation and variance in a coupled pair of processes started from fi^. 

Proposition 2.1. Let w(-) andtu^(-) be a coupled pair of TAZRP started from 
the distribution fi^. Let Q a (-) be a second class antiparticle between the processes 
starting at Q a (0) = 0. Let the model evolve according to the basic coupling. 
Then the variance of the current for any V G M in the stationary process can be 
rewritten as: 

VarV°°(*)) - Var<>ME( | [Vt] Q a (t)\) 
= g(l + g)E(\[Vt]-Q a (t)\). 
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Also, 

(2-6) E(Q a (t)) = VH=j^-^t 

holds. 

In article [5] identities (2.5) and (2.6) are proved for a class of attractive 
models that include the asymmetric simple exclusion, zero-range and bricklayer 
processes. 

The interesting current fluctuations occur at the characteristic speed V e = 
(i+ e )2 • From (2.6) we see that this is the average speed of the second class particle. 
The characteristic speed also appears in the PDE theory of the conservation law 

(2.7) Qt + f(g) x = 

that is obtained via the Eulerian hydrodynamic limit of TAZRP. Here the hydro- 
dynamic flux takes the form 

(2-8) 1(g) - E«(I (a , >0) ) = T i-. 

At constant density g the characteristic speed is f'(g) — V e at which perturba- 
tions of entropy solutions of (2.7) travel. For more information we refer to [10] 
for hydrodynamic limits of interacting particle systems. 

Let us return to the stationary TAZRP with Geometric (jj^) occupation vari- 
ables at each fixed time. A basic object for understanding space-time correlations 
is the two point function S (i , t) = E e (oji(t)ojo(0)) — g 2 . Due to [5] it can be written 
as the transition probability of the second class particle: 



(2-9) 




= Var e {cj a )-P{Q(t) = n}. 



The sum under the expectation is regarded as a Radon-Nikodym derivative that 
takes our second class particle in the initial setting p, where the second class 
particle resides initially at the origin. Note that use of the term "transition 
probability" is not meant to suggest that Q(t) is a Markov process. 
From (2.9) and (2.6) follows 

*S(i, t) = Var e (uj )V Q t = g(l + g)VH. 

The diSusivity is by definition a normalized second moment of the two-point 
function: 

S« = ^I>-V'<) 2 S(M>. 

Consequently the diffusivity can also be expressed in terms of the variance of the 
second class particle: 

D(t) = t- 1 Var(Q(t) ). 
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Below we show that, similarly to the ASEP, the TAZRP second class particle is 
also supcrdiffusive with variance of order i 4 / 3 . 

We can now state the main theorem, a moment bound for the second class 
particle in a pair started from fi^. 

Theorem 2.2. For the constant rate TAZRP started in fi^ distribution for 
any g > there exist constants < to, C < oo such that for real 1 < m < 3 and 
t > t Q 



(2.10) C- 1 < E 



Qa{t)-VH 



t 2/3 



< C. 



Let us simplify notation to J e (t) — J^ v \t) for the current across the char- 
acteristic. Via (2.5), taking m — 1 gives the variance of this current. 

Corollary 2.3. With assumptions as in Theorem 2.2, for t>to the variance of 
the current across the characteristic satisfies 

C -i t 2/3 < V ar e (J e M) < Ct 2/3 . 
Taking m = 2 identifies the order of the diffusivity. 
Corollary 2.4. With assumptions as in Theorem 2.2, for t > t 

C-h 1/3 < D(t) < ct 1/3 . 

The upper bound of (2.10) is not valid for all t > 0. For small t the variance 
of Q(t) is of order t because the second class particle is likely to have experienced 
at most one jump. 

The above results are completely analogous to those in [4]. 

Next, another important consequence of the theorem is the Weak Law of Large 
Numbers for Q(t). Markov's inequality implies 

Q(t) ye _^E|Q(t)-[^]| + l^,CT 2 / 3 



t 



> e ) < 1 u < > 0. 

et et 



We mention here that the moment bounds of (2.10) together with the fact that 
Q(t) is dominated by a Poisson process are strong enough to derive the Strong 
Law of Large Numbers as well. 

Proposition 2.1 and Theorem 2.2 imply that the current variance is of order 
* for all V ^ V e : 

Var e jW(i) i 
lim 1 \H = lim -Var e (w )E(|g(t) - [Vt]\) 

t-^oo t t—+oo t 



= lim ivar e M(E|Q(t) - [VH] + [VH] - [Vt]\] 

t— >oo t \ J 

= Var e (u> )|^ e - V\. 
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The Central Limit Theorem for J^ v \t) also holds in the following form: 



/ J ( y )(t) WJ^jt) \ $ , x for aU R V^T/e, 
V*Var e (w )|ye-^| / w 

where $(x) denotes the standard normal distribution. 

A more important consequence of Theorem 2.2 concerns the fluctuations of 
the tagged particle of the TASEP near the characteristics. To be able to ex- 
plain the result we introduce some new notation: first, wc consider a stationary 
TASEP with independent Bernoulli (a) occupation variables at each site which 
we condition on the event that the origin is occupied. Next we label the parti- 
cles in an increasing order such a way that the one at the origin becomes label 
0. We denote the site of the i-th particle at time t by Rf(t). Then we let the 
system evolve according to the TASEP evolution with particles jumping to the 
left and not to the right. Rates of these jumps are 1{ u e =1 } • I^e _ j, i.e. the 
particles are allowed to jump one to the left if the neighboring site on the left is 
empty. The particles in the TASEP cannot pass each other, hence the ordering is 
preserved over time. As we conditioned on the origin being occupied, i2jf(0) = 
and we choose the usual height representation of the process with Hq(0) = and 
the height sloping to the left. Since the particles jump to the left, this height 
function is decreasing in time. 

Second, we look at the process of the differences, namely the empty sites 
between two neighboring particles in the exclusion process. Define u>i(t) := 
Rf (t) — Rf_ l (t) — 1 . Then the process w (t) is a stationary TAZRP with parameter 
q = l/a — 1, where particles jump to the right. The height representation hit) 
of the zero range process can be translated to the TASEP, namely 

(2.11) R%(t) = -h k (t) + k. 

For k = [VH], wc get Rf VH] {t) = -J e (t) + [V e t\. Theorem 2.2 leads to 

Var(i?^ (t)) Var(fljs (t)) 
C- 1 < liminf -4tt^ — - < limsup -^rr 1 — ~ < C 

To give an interpretation of this result, we recall that the characteristic speed 
in the TAZRP equals V s = nj^p = " 2 ' so equation (2.11) at site [V e t] turns 
into R [a 2 t] (t) = -J e {t) + [a 2 t]. Since the mean distance between particles in 
the TASEP is l/a, the [a 2 i] th particle started at time somewhere near at. 
Burke's Theorem implies that each particle marginally moves to the left following 
a Poisson process with parameter 1 — a, hence this particle should be at time 
t somewhere near (2a — l)t, which is exactly the characteristic position of the 
exclusion process. So, to put the result Var(i?^ 2t j) ~ t 2 ^ in words we can say 
that the position of the tagged particle that is expected to be on the characteristics 
has fluctuation of order t 1 / 3 . 

The rest of the paper is devoted to the proof of Theorem 2.2, first the upper 
bound and then the lower bound. 
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3 Upper bound 

Proof of the upper bound utilizes couplings of several processes. Before deriving 
the upper bound we introduce some preliminaries on the couplings. 

Couplings 

We start by describing an initial distribution of three coupled processes and with 
labels attached to second class particles between two of the three processes. 

Definition 3.1. For two probability measures < fi 2 denotes stochastic domi- 
nation. 

Equivalent to stochastic domination are: 

• For the cumulative distribution functions F\ (x) > F 2 (x) holds for all 

• There exists a coupling measure fj,(x, y) with marginals \i\ and \i 2 such that 
fj,{{x,y) : x < y} = 1. 

Fix densities < A < g. In the next sections we want to couple models with 
different densities and distribution so here we prove stochastic domination of the 
different measures. 

Lemma 3.2. The measures [i e and fi e are monotone in g, i.e. 
1- Ma < Me f or a M A < f?, 

2. fix < fi e for all A < g 

3. Furthermore, fi e > /U e for all g > 0. 

Proof. The cumulative distribution function of the geometric distribution with 
mean g is 

Differentiation with respect to g gives that the function is monotone decreasing in 
g, meaning that A < g implies F e (z) < F\(z). Next, fi g is the distribution of the 
sum of two independent, geometrically distributed random variables, implying 
Fg{z) < F\{z). The third statement comes from the fact that a random variable 
of distribution fi e is the sum of two independent random variables of distribution 
Me- □ 

Definition 3.3. 

• v is the coupling measure with marginals [i\ and [i e with v{{n,uj) : rj < 
oj} = 1. Let v_ := v® 1 be the product measure with marginals v for each 
site. 
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• is the coupling measure with marginals fi\ and fi e with 0{(r],u)) : 77 < 
w} = l. 

We now describe the initial distribution. Let (77, w~) be distributed in the 
product measure of marginals v for each site i ^ and v at the origin. We also 
couple a third process denoted by w to these two processes with ui(0) = oj~ (0) + Sq 
i.e. we add an extra second class particle to the process w~ . Saying differently we 
have a second class antiparticle on u. Let v_ denote the resulting joint distribution 
°f {Vi <^~> <±0- Then £-a.s. we have the following initial conditions for the processes 
77, oj~ and uj : 

• V < w~ < u, 

• There is a second class particle between u and w _ at the origin at time 0, 
implying that there is at least one second class particle between uj_ and 77 at 
the origin. 

• There are infinitely many oj — 77 second class particles on both sides of the 
origin. 

Label the lo — 77 second class particles with integers in an increasing fashion 
from left to right, giving label to the topmost second class particle initially at 
the origin. Xk(0) denotes the position of second class particle with label k, so 
that 

(3.1) • • • < X_!(0) < X (Q) = < Xi(0) < • • • 

Let these configurations evolve from the initial distribution v_ with common 
Poisson clocks. The rule for the labeling is that whenever a second class particle 
jumps from site i, it is always the highest labeled one at i and it arrives at the 
lowest place at site i + This way the ordering (3.1) is kept for all time, and 
the basic coupling preserves the ordering rj(t) < u~(t) < uAt) a.s. for all times 
t. The topmost u — 77 second class particle that starts at the origin is of special 
importance to us, so we set X{t) = Xo(t). The effect of the coupling is that the 
u> — lu~ particle (denote its position from now on by Q a {t)) cannot pass X(t). 
Precisely saying, the ordering Q a (t) < X(t) is preserved over time. The argument 
for this fact is that initially Q a {0) < X(0) holds. Since both particles are only 
allowed to jump to the right Q a can pass X if and only if they are at the same 
site. In this case X(t) jumps if rjx(t)(t) = but Q a {t) jumps only if ^ = 0, 
and the coupling 77 < lu~ guarantees that the second event is part of the first one, 
i.e. the jump of Q a (t) implies the jump of X(t) but not conversely. 

Throughout this section probability P refers to the law of the three processes 
started in distribution z>, evolving according to the basic coupling as described 
above. 

Proof of the upper bound 

We begin with the proof of the upper bound in Theorem 2.2. C and its variants 
C\ , C2 , ■ ■ ■ denote positive constants that possibly depend on g and whose values 
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can change from line to line. We first prove a lower bound on the density of the 
uo — rj second class particles. For integers jeZ and u > let 

j'+2«-l 

(3.2) N j (t)= J2 M*) "%(*))■ 

i=j+i 

Lemma 3.4. Let g > A > and d > be an integer. Then there are strictly 
positive finite constants 7 = 7(g), Ci = C\{g,d) and C'2 = Cz(q) such that the 
following holds: if < g — A < 7, t/ien for all integers j e Z, it > and any time 
i > 0, 

p{Nj(t) < u{g - A) + d} < Ci exp{-C 2 u((? - A) 2 } 

Proof. For the moment, denote by w(-) a constant-rate TAZRP such that yi(0) = 
u>i(0) at all sites i except for i — 0, and z(-) is a constant-rate TAZRP process 
such that Zi(0) = r/i(0) at all sites i except for z = 0. For i = 0, we pick the 
pair (zo(0), yo(0)) in distribution v, independently of the configuration on other 
sites, such that z(0) < y(0) < w~(0) and z(0) < n(0) < w~(0) holds. This can 
be done since /ia < Me < fi e and /ia < Aa < As- Apply the basic coupling to 
ensure < 77(f) < w(f) and y(f) < w~(f) for all f > (notice that this holds 
initially). y(t) and z(t) are marginally time-stationary processes, hence we can 
omit the notation for their time dependence in our arguments. However, the pair 
(z(t), y{t)) is not in product distribution for t > 0. Define 

j+2u-l j+2tt-l 

F = ^ j/i and Z = ^ 

i=j+l i=j + l 

so that JV^t) = EiZjZi\"i(t) -y*{t)) + Y-Z + E^li'iMt) - rh(t)). In this 
expression the two sums can be estimated from above by the total number of lo — y 
and z — 77 second class particles, respectively. Using the facts that these models 
only differ initially at the origin and the total number of second class particles 
is preserved over time by the basic coupling allows us to use the estimations 
< £^i~W)-y*W) and Eg^teW "WW) > *o(0)-»jo(0) > - % (0). 
If we plug in the lower bounds the probability of the event increases: 

P{jv,-(t) < u(g - A) + dj < P{r -Z- ffo(0) < u (e — A) + d| 
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Now we use exponential Markov's inequality to get for any a > 
p[Nj(t) < u{q- A) +dj 

< e au(g-\)+da _ Q aZ e ar)o(0)^ 

< e ««(<?-A)+d Q . [ E <?(e- 4 « y )] 1/4 • [E A (c 4aZ )] 1/4 • [E(e 2a, »< >)] 1/2 
_ e c«((?-A)+<i Q . [Ee(e _4a! " ) )] 1/2u_1/4 • [E A (e 4azo )] 1/2 " _1/4 

. [E^o (0)^1/2 

= exp{au(e - A) + da - (^u - J) log[l + g(l - e~ 4tt )] 

- (1« - 1) log[l + A(l - e 4 ")] - log[l + A(l - e 2 ")] } 

< exp{-au(g - A) + 4u(g + A + £ 2 + A 2 )a 2 
+ MC 3 a 3 + C 4 a + C 5 a 2 + C(a 3 )}. 

Here we used the marginal fj,^ and product distributions of z and y, and 

the fact that 770 (0) ~ A a which implies E(e 2m)o(0) ) = (1 + A + Ac 2 ")" 2 . The last 
inequality comes from Taylor expansion w.r.t. a. The 0(a 3 ) term is uniform over 
A < g for a fixed £>. In order to ensure the existence of the exponential moment we 
need the condition j^e 2a < 1 to be satisfied, or equivalently, a < i log(-Mp-). 
Next we pick 

g-A 

° 8(e + A + £ 2 + A 2 )' 

optimizing the g and A-dependent terms in which u is also present. Comparing 
this to the required condition above gives a bound for 7, i.e. 7 < 4 log ( — ^) ■ (f? + 
p 2 ). By this choice of a we get 

P{^(t) < «(, - A) + d} < e«^^<» S » exp { - 16(g ^;y +A2) } 

Now g — A < 7 gives rise to the constant Ci in the first factor and finishes the 
proof. □ 

Now we are ready to turn to the main estimate for the upper bound. The idea 
is to bound the deviation P{Q a (t) > u + [V e i]} with an appropriate expression 
that involves the moment E\Q a (t) — This is completed in (3.7) below and 

then the upper bound comes from an elementary integration step. 

Along the way we compare currents in two processes that we abbreviate as 
follows: 

J e {t) — J x +e A ' (t) for current in the process, and 
J ye ' x (t) = t+i?\t) for current in the r}(.) process. 
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Notice that both use the same speed V e = f'(g) = j^t for the observer. As 
already defined in the Introduction, this is the characteristic speed of the 
process with density g. From now on we denote by tilde the centered random 
variable. 

Lemma 3.5. Suppose g — A < 7 with 7 from Lemma 3.4- Then for positive 
integers u and times t e [0, 00), 



P{Q a (t)>2u+[VH}} 



2 



(3.3) < P{J?(t) - J VQ '\t) > u(g - A) - } 

+ dcxp{- C 2 u(g- A) 2 }. 

Proof. Now we recall the construction of the labeled particles and the fact that 
Qa(t)<X(t) Vt: 

P{Q a (t) >2u + [VH]} < P{X(t) >2u + [V e t]}. 

N[vet]{t) 01 (3-2) counts the number of ui — r\ second class particles at time t in 
the interval {[V^ e £] + 1, . . . , [V t] + 2u— 1}. Since the second class particles stay 
ordered and X{t) started at the origin, the event {X(t) > 2u+ [V e t]} implies that 
all these second class particles crossed the path s 1— » [V e s] + 1/2 by time t. Each 
such second class particle crossing increases J s (t) — J ye ' X (t) by one. Therefore 

P{X(t) >2u+ [VH]} < P{J e (t) - J ve > x (t) > N [vet] (t)} 
< P[j e (t)-J ve ' x (t) > u(g-\) + 4g+l} + p{N [Vet] (t) < u(g- A) + 4g+ l}. 

Combine the previous displays to get 

(3.4) P{Q a (t) >2u+ [VH]} < P{.J e (t) - J ve > x {t) > u(g - A) + Ag + l} 

(3.5) +p{N [vet] (t) <u(g-\)+4g+l} 

To line (3.5) apply Lemma 3.4 with d = \4g + 1] . We see that line (3.5) is 
bounded by the last exponential term in (3.3). Recall the definition of the flux 
(2.8). If lo and 77 start from their respective n and n x equilibria then we would 
have, due to (2.4) and (2.6), 

E«(J C (*)) - E x (J ve - x (t)) = t(f{g) /(A) - f'(g)(g A)) 

(l + g) 2 (l + A) 

where we ignored the error coming from the integer part of V e t. This error can be 
estimated from above by the number of particles sitting at one site in each of the 
processes, i.e. W{je{t)) - E x (J vQ,x (t)) - 1 ( 4$(i+X) ^ ^ e (^) + ^(m) = Q + >>- 
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Our processes are also perturbed initially at the origin by v being different of v, 
which gives an error E(J e (t)) - E e (J g (t)) < E(w o (0)) = 2g+l. In the 77 process 
the term —E(J v ^ x (t)) + E x (J v e' X (t)) is negative so it can be estimated from 
above by 0. The term Ag + 1 inside the probability on line (3.4) makes up for 

these errors. So, we get P{ je(t) - J v "' x (i) > u(g - A) - t (1+ ( ^ A j^ A) } Finally, 
we omit the factor 1 + A in the denominator, implying the decrease of the right 
hand side of the inequality, and increasing the probability. □ 

Lemma 3.6. 

Var(J e (t)) < 120(0+1) + 2 + 2g(l + g)E(\[V s t]-Q a (t)\), 
Var(J v ' A (t)) < 12g(g + 1) + At(g - A) + 2g(l + g)E(\[VH] - Q a (t)\) . 

Proof. The variance Var in the statement is taken in the three-process coupling 
where the geometric distribution fj, is initially perturbed at the origin. Recall 
that Var e denotes variance in the stationary process y(t) coupled to w(t), and 
denote by J|(i) the height function in this process at site [V e t] at time t. The 
following estimation holds: 

V&r(J e (t)) < 2Var(J e (t) - Jj>(t)) + 2Var e J e (t) 
(3.6) < 2E(co (0) 2 ) + 2Var e (^ (0)) • E(|Q B (t) - [VH]\) 

< 12g(g+l) + 2 + 2g(g+l)E(\Q a (t)-[V e t}\). 

Here we applied Proposition 2.1 to the second term Var e (J e (t)) in the first line 
and used the fact that J e {t) and Jy(t) only differ at the origin by Wo(0) — yo(O), 
so its variance can be estimated from above by the second moment of Wo(0), 
equal to the second moment of 1 plus the sum of two independent geometrically 
distributed random variables with mean g. Similar computation for the A-density 
process is as follows (with stationary A-density process z(t) and its height function 
JY e '\t)): 

Var(J v ^(t)) < 2Var(J v '^(t) - J z y ^ A (t)) + 2Var A J v ^ x (t) 

< 2E( % (0) 2 ) + 2Var A (r7 O (0)) • E(\Q x (t) - [VH]\) 

< 4A(A + 1) + 8A 2 + 2A(A + l)E(|Q A (t) - [VH]\) 

< I2g(g + 1) + 2A(A + l)E(\Q x (t) - Q a (t)\) 

+ 2g(g+l)E(\Q a (t)-[VH}\) 

< I2g(g + 1) + A{g - \)t + 2g(g + l)E(|Q„(t) - [VH]\). 

In the first line we cut the variance into two terms, where the first term can 
be estimated from above by 2E(?7o(0) 2 ), and for the second term we applied 
Proposition 2.1 with a second class particle Q x (t) added to the 77 process. Then 
we used triangle inequality and Proposition 2.1 again to get E(\Q x (t) — Q a {t)\) = 
E(Q x (t) — Q a (t j) = t mj^p ~ ^ 7TTeP ' a b som te value can be omitted by 

Q X {t) > Qa(t), which we show below. Initially Q A (0) = Q o (0) = 0. The jump 
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rates of Q x (t) and Q a (t) are \ Q \ {t) (t)=o an d I w - (t)=o> res P ec tively. Second class 
particles are only allowed to jump to the right, and a right jump of Q a without 
Q x is impossible by the basic coupling and the above rates since r/j (t) < uu^ (t) . 

Some algebra then leads to the estimation 2A(A + l)E(|Q A (t) — Q a {t)\) < 
A{g — \)t. Collecting terms completes the proof of the lemma. □ 

We come to the lemma that summarizes all the previous estimations. 

Lemma 3.7. For any real u > 1 and time t > 0, 

(3.7) P{Q a (t) >2u + [VH]} < C 3 ^ • E(|Q„(t) - [VH]\) 

+ c V + c V + Ciexp {- C2 7J} + e_tl - 

Proof. Set b = jj^p (lAg) where 7 is the constant from Lemma 3.4. We proceed 
by considering three cases. 

Case 1: 1 < u < bt. Suppose first u is an integer as it was in the proof of 
(3.3). Throughout density g has been fixed, and now we also fix 

The constraint on u guarantees that A > and g — A < 7 which was required 
for (3.3). The point of this choice of A is to maximize the lower bound inside 
the probability on the right-hand side of (3.3). So, continuing from (3.3) with 
Chebyshev's inequality and the previous lemma, 

P{Qa(t) > 2u+[VH]} 

< Pp(t) - J^ x (t) > £(1 + g) 2 } + C lC xp{-C 2 ^} = 

< C 3 ^Var(J^) - J ye ' X (t)) + Ci cxp{-C 2 ^} 

< C 3 ^(2Var(J^) + 2Var J V °> X {t)) + dcxpj-C^} 

t 2 t 2 t 2 f ?; 3 1 

< C ^ ■ E(|Q-(t) - W e t]\) + G,- % + C 5 ^ + & exp{-C 2 ^}. 

Extension from integral u to real u is achieved by adjusting constants on the last 
line above. 

Case 2: bt < u < It. Then bu/2 < bt and 

P{Q a (t) >2u + [VH]} < P{Q a (t) > 2 • bu/2 + [VH]}. 

Case 1 can be applied with u replaced by bu/2, at the price of adjusting some 
constants with powers of b. 
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Case 3: u > 2t. Since Q a (t) is bounded above by a rate one Poisson process 
N(t), 

P{Q a (t) >2u + [VH]} < P{Q a (t) > 2u} < P{N(t) > 2u} 
< e- 2u • E(e JV < t )) = e t(e - 1 '>- 2u < e~ u 

for all t. Combining the bounds from the cases proves the lemma. □ 

We are ready to complete the proof of the upper bound of Theorem 2.2. 
Introduce the notation 

*(t)=E(|Q (t)-[V«i]|). 

In order to get a bound on P{\Q a (t) — [V e t]\ > 2u] we need a lower tail 
bound for Q a (t), similar to (3.7). This can be obtained by methods analogous to 
the ones we have applied throughout Section 3. In this case we have to couple a 
process with initial distribution fx for some a > g close enough to g and derive 
similar lemmas to get the same conclusion for P{Q a (i) < — 2u + [V e t]} as in 
(3.7) with readjusted constants. 

Introduce a large constant 2 < r < oo. For t > 1 and u > rt 2 ^ 3 we can 
combine (3.7) and the matching lower tail bound. Replace 2u by u (wc made 
sure u/2> 1), then some algebra leads to 

P{\Q a (t)-[VH]\>u} 
(3 ' 8) < + C 2 (r) ( J + exp{-C 3 (r)^3 } ) . 

The two exponential terms in (3.7) were combined via e~ u < exp(— ut~ 2 / 3 ) and 
exp(-C 2 u 3 t- 2 ) < exp(-C 2 r 2 wt- 2 / 3 ). 

Let 1 < m < 3. Integrate bound (3.8) over u e [rt 2 / 3 , oo): 

E(\Q a (t)-[VH]\ m ) 

OO 

(3.9) <r m t 2m / 3 +m J P{ \Q a (t) - [VH] \ > uju^ 1 du 

rt 2/3 

<C 1 r m - 4 f(f)t 2+ ( 2 / 3 )( m - 4 ' + C±(r)t 2m l 3 . 

The constant C\ depends on m but not on r, while C4 depends on both m and 
r. To get the final bounds, take first m = 1 in (3.9) to get 

*(t) < C ir - 3 ^(t) + C 4 (r)t 2 / 3 . 

Since C\ is independent of r, fixing r large enough gives ^(t) < C^{r)t 2 / 3 . Using 
this bound in the last line of (3.9) implies 

E(\Q a (t)-[V s t]\ m ) <Cs(r)t 2m ' 3 

for 1 < m < 3, which completes the proof of the upper bound of Theorem 2.2. 
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4 Lower bound 



The lower bound is proved by perturbing the stationary distribution of a process 
on a segment of the lattice. We start with discussing the initial distribution of 
some coupled processes. 

Perturbing a segment initially 

Recall again the characteristic speeds V e = nqr^p and V x = njxp ■ We set 
q > A, hence V e < V x . Throughout this section u > denotes a fixed positive 
integer, and 

n = [V x t] - [VH] + u. 

Recall the Definition 3.3 of the measures v and v. Define an initial product 
distribution of two configurations (77(0), C(0)) by setting the marginals on each 
lattice site i: 

{ (T?i(0), Ci(0)) ~ v Xi<-n, 

(77,(0), Ci(0))~i> if* = -n, 

rn(0) = Ci(0) ~ MA if-n<«<0, 

LM0),Ci(0))~i/ if * > 0. 

Note that the number of particles at site — n is fi\- distributed for r;(0). Except 
for this perturbation 77(0) starts in the stationary Geometric product distribution 
fly The process ((■) initially has distribution fi , except at sites {— n+ 1, . . . , 0} 
where the parameter g has been replaced by A, and at site — n where it has 
measure /t e . 

We add a second class particle to the process rj(.) initially started at site — n 
and denote its position at time t by Q(~ n \t). Introduce, as before, ri + {t) := 
V(t)+S Q -„ {t) . 

Define a third initial configuration by 



6(0) 



Ci(0) ifi<-n, 
J?i(0) if i > -n. 



Apply the basic coupling to obtain the joint evolution of all these processes. This 
guarantees the majorizations 

nit) < l(t) < C(t) and h c (t) < h*(t) 

where the last inequality is for column heights. 

As before in (2.3), denote the net particle currents by J V ' n and J v ^ in the 
respective processes jj(-) and ((■). The first observation is that Q^~ n "> gives one- 
sided control over the difference of these currents. 

Lemma 4.1. For any V € M 

Q(-")(t) < [Vt] impKai J y ' c (t) - J v '"(i) < 0. 
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Proof. Denote the positions of the ^ — 77 second class particles at time t by 

■ • • < Y k (t) < ■ ■ ■ < Y_ 2 (t) < < Y (t). 

This order of the labels is again preserved over time. Initially Y o (0) < —n = 
Q(-«)(0), and the respective jump rates I VYqW =o and Ir )Q( _„ )(t) =o of the second 

class particles ensure Y)(t) < Q^~ n ^(t) for all later times. Indeed, from the time 
Y has jumped on Q^~ n ^ they always jump together. So, Y cannot leave Q^~ n ^ 
behind itself. Saying it another way, once Q(~ n \t) = Y a (t), this property will 
hold forever. 

Now note that there are no £ — r] second class particles strictly to the right of 
Q(~ n \t) at time t. So, the height difference is zero to the right of Q^~ n \t). Also 
recall the majorization ft^(i) < ht(t). Thus under {Q^^it) < [Vt]} we have 



= hf vt] (t) - y vt] (t) > rf vt] (t) h? vt] (t) = J v -Ht) - J v «(t) 



"[Vt] 



□ 



Let w(-) be a TAZRP started from the distribution ft shifted n sites to the 
left, i.e. at site — n the number of particles is /^-distributed. The next lemma 
compares the distributions of ( and u. 

Lemma 4.2. Denote by P w and the probability of events that depend only 
on the respective processes (!>(•) and £(•). Then there exist 7 > such that for all 
A > g — 7 the following inequality holds: 



pC(.) < P*(.)5 . ex p 



n(g — A) 2 



g(l + g) 



Proof. Let 



z= £ Ci(o). 

i— — n+1 

Z has a Negative Binomial(n, j^) distribution, namely 

^™-*>-(T70(m)'(rh)" 



for z > 0. We use the Cauchy-Schwarz inequality below to perform a change of 
measure on this negative binomial distribution. 



z=0 



00 



' . > JP c (-|Z-z)] 2 m e (z) 

~z=0 



m x (z) 
[m«(z)]5 

\ OO 

'■it;- 

z=0 



m,e(z) 
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Now the last part equals £~o (^V) " ( 

can recognize the negative binomial mass function with parameter 
< 1 — Ci+TFg < 1, so we use the identity 



In this formula we 



E 



-l\/ A 2 (l + g) N* / 

1 \(1 + X) 2 q) \ 



(1 + A)V 



A 2 (l + g) 
(1 + A) 2 g 



= 1 



and the remaining factors equal 
So we get 



A 2 (e+i) 



• ( s+L. 



»(•)< [f><(-|Z = *)] 2 m«(*) 

n i 

< \£vt(.\Z = z)m°{z) 



1 - 



(g-A) 2 



2=0 



• exp 



0(1 + 0) 

n(g-A) 2 



g{\ + g) 



In the last line we used the fact that l — x>e 2x for all < x < 0.5. This will be 
satisfied for x = , for A close enough to g. All that is left is to recognize that 

P^( • | Z = z) is the probability depending on a process £(•) whose initial distri- 
bution coincides with the distribution of Cj outside {— n+ 1 . . . 0}, with z particles 
distributed in that interval with each configuration equally likely. Namely, each 

configuration with z particles in this interval has probability (t+x) (i+~a) ' so 
conditioned on the event that z particles are here initially, each configuration 
has probability ■ Summing these conditionals, weighted with the Nega- 

tive Binomial(n, yq-;) mass function m e (z), gives the product Geometric initial 



distribution n of the process £>(•). 



□ 



Proof of the lower bound 

Now we only need two steps to obtain the lower bound. In this part we need 
again the concept of a second class antiparticle started from the origin on a \i - 
stationary process initially perturbed by setting u>o(0) ~ jl g + So- We denote 
its position at time t by Q a {t). The quantity of primary interest is abbreviated, 
as before, by = E(\Q a (t) - [VH]\). 

The first step is proving an upper bound on Q^~ n \t). As before u is an 
arbitrary but fixed positive integer and n = [V x t] — [V e t] + u. 

Lemma 4.3. 

(4.1) P{Q^(t)>[VH]}<^ + 4 -^^ + 2 -. 

u u u 
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Proof. Below Q x (t) stands for the position of a second class particle started from 
the origin on a process rj(-) in fi^ distribution. Translation invariance implies 

P{Q<" n >(t) > [V e t]} 

= P{Q ( - n) (t)+n-[V x t] >u} 
= P{Q x (t)-[V x t]>u} 
^ E(\Q x (t)-[V x t}\) 
u 

< E(\Q x (t) - Q a (t)\) + E(\Q a (t) - [VH}\) + [VH]~[Vn} _ 
u u u 

As in Lemma 3.6, the first term equals ^ ( (j+yp — bounded from above 

by -t(g — A) after some calculations. The second term is ^(t)/u, and the third 
term is similarly estimated by -t(g — A) + -, the last part coming from possible 
integer part errors. □ 

The second step is an estimate of the probability of the complement of the 
event in (4.1). 

Lemma 4.4. For any < K < t ^'^l - 2g, 

2g(l + g) (6 + *(*)) 4t(g-\) 
K 2 K 2 ■ 



Proof. Lemma 4.1 leads to 

P{Q<- n >(t) < [VH]} < P{J v ^(t) - J v ^(t) < 0} 

(4.2) < P{ J™(t) <K + t( T A_ _ _^_) + 2A} 

(4.3) + F {jV^ t)>K + t{ ^__^ )+2x} 
We apply Lemma 4.2 to line (4.2) to bound it by the probability of the process 
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(4.2) < [P"{J e (t) < K + t(- 



X 



X 



< [P ffl {.7«(t) <K-t— 



+ X (1 + qY 

{e-xf 



) + 2A}] 5 • exp 



n(g-X) 2 



(l + X)(l + g) 2 

<[P*{j«®<K-tte^ + 2e}Y 

n(g-X) 2 



+ 2A}] 



exp 



g(l + g) 
n(g - A) 2 " 



exp 



n(g-A) 2 



< Kgffl^! -exp 



+ (e-A) 2 
TT+eF 



2g-K 

1/2. 



< 



M 2 



L (e-A) 



- 2q- K 



t (e-A) 2 
(i+e) 3 



/a "I 

-2g-K> 



exp 



n(e- A) 2 



g(l + g) 



When centering, we used J e (t) = J e {t) + [EJ e (t) - E e J e (t)} + [E e je(t) -t(j^- 
tt+op )] ~*~^(l+g ~ (i+g) 2 ) ' wnerc both error terms in the brackets [•] are nonnega- 



tive. So, JB(t) <K + t(j^- ( I ^)+2Aimplies J«{t) < 



(g-A) 2 



i+A (i+ e ) 2 / 1 -'"—K"^" vv - -"■ "(i+A)(i+ e ) 2 r2A - 
Then, in the last line we used Lemma 3.6 to bound the variance by the function 
ty(t) even though the second class particle in u) starts at — n rather than at the 
origin, namely Var(J e (i)) < 2 + 12g(l + g) + 2g(l + g)^!(t). We use the inequality 
\Ja + b < v / a + Vb in order to get rid of the term +2. The only change needed in 
the proof of Lemma 3.6 is in the calculation (3.6) where we must add E(u_„(0) 2 ) 
instead of E(o> (0) 2 ), but these are equal. 

A similar coupling consideration shows that E(J y ^ r '(t)) differs by at most A 
from the same expectation taken under a stationary /x A initial condition. Thus 
taking integer parts again into account, giving another error term A, line (4.3) is 
bounded from above by 



P{j ve ' r >(t) > if} < 



Var(J v "■") 
K 2 



Lemma 3.6 can be applied again to bound this variance with the similar change 
of shifting all processes by n sites to the left. Hence we can continue from above 
to bound line (4.3) with 

2<?(1 + e) (6 + *(*)) 4*0? -A) 



K 2 



K 2 



□ 



Now we are at the last step of proving the lower bound of Theorem 2.2. By 
Jensen's inequality it suffices to prove for the case m = 1, in other words that 



liminf t- 2 'H{t) > 0. 

t— >OG 
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In the last two lemmas take 



u = \ht 2 ^] , q - A = fa" 1 / 3 , and K = bt 1 / 3 , 

1 l2 1 

where h and 6 are large, in particular b large enough to have bt 3 < ( 1+g )a ^ — 2g 
so that if satisfies the assumption of Lemma 4.4. Then 

n = [VH] - [VH] +u< —l±^±^—{ e -X)t + 2 + u<(2b+ h)t 2 ' 3 < Ct 2 ' 3 

for large enough t. We can simplify the outcomes of Lemma 4.3 and Lemma 4.4 
to the inequalities 

(4.4) P{Q ( - n »(i)>[B]}<cf 4 + 



£2/3 h ht 2/3 

and 

1/2 

t l/3- 



The new constant C depends on b and h. 

The lower bound now follows because the left-hand sides of (4.4)-(4.5) add up 
to one for each fixed t, while we can fix b large enough and then h large enough 
so that 4b/h + 4/6 < 1. Then t^ 2 / z ^{t) must have a positive lower bound for all 
large enough t. This completes the proof of Theorem 2.2. 
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